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INTRODUCTION 

Hadronic collisions in the high energy limit are probing with a better and better time resolution the quantum fluctua- 
tions in the wave-functions of the incoming particles, due to time dilation. On the contrary, due to length contraction, 
the longitudinal resolution is getting worse and worse. In the Regge limit, the typical transverse scales related to a 
given scattering process are kept fixed while the high energy limit of the collision is taken, so that the transverse 
resolution stays constant. Hence, in the Regge limit, a growing number of fluctuations of the incoming particles are 
participating to scattering processes, but those located in an volume set by the transverse and longitudinal resolutions 
cannot be distinguished and act coherently. Coherent multiple scattering of quarks and gluons is thus an unavoidable 
ingredient in the physics of the Regge limit. Equivalently, the coherent clusters in the hadronic wave-functions are 
more efficiently described in terms of semi-classical gluon fields than of partonic Fock states. 

The physics of the Regge limit is then beyond the reach of the improved parton model, with parton distributions 
satisfying the DGLAP and/or BFKL evolution equations, in which the hadronic wave-functions are assumed to be 
dilute, and scatteringprocesses are assumed to occur via a single binary collision of partons. The phenomenon of 
gluon saturation |[lll3 corresponds to that breakdown of the improved parton model picture via the onset of high 
density effects. Gluon saturation as well as the following dense regime are properly described in the Color Glass 
Condensate effective theory [3, 4, 5] and related formalisms. 

These frameworks have been successfully applied to the calculation of various observables, related for example to 
forward physics at RHIC and at the LHC, and to early stages of heavy ion collisions. But the most straightforward and 
best understood application remains the one to deep inelastic scattering (DIS) at low Bjorken x. There is an ongoing 
effort for a few years to update from leading order (LO) and leading logs (LL) to next-to-leading order (NLO) and 
next-to-leading logs (NLL) the theory of DIS with gluon saturation. The purpose of the present contribution is to 
review the main progresses in that direction. But first, let us recall the LO/LL formalism and results. 



GLUON SATURATION AT LO/LL ACCURACY 

The Ft and DIS structure functions are proportional to the total cross sections for the collision between the target 
and a transverse or longitudinal virtual photon. At low x and at LO accuracy, these total cross sections satisfy the 
dipole factorization fS): they can be written as a convolution of the probability of splitting of the photon into a quark- 
antiquark dipole before reaching the target, known from QED, with the imaginary part of the dipole-target elastic 
scattering amplitude, which is a non-perturbative quantity. 
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However, logarithmically divergent soft gluons emissions by the dipole are part of the higher order corrections. 
The corresponding leading logs can be resummed via a renormalization group evolution of the dipole-target scattering 
amplitude with respect to the total energy of the collision or to the available rapidity range Y ^ log( 1/x). This evolution 
in Y is given by the B-JIMWLK equations |7, 8, 9, 10, 11, 12, 13, 14, 15], which, for the DIS case, can be safely 
approximated by the BK equation ITI I 1611 1711 . which is a nonlinear extension of the BFKL equation. 

The nonlinearity of the BK equation encodes gluon saturation and restores the unitarity of the dipole-target scattering 
amplitude at fixed impact parameter. The impact parameter dependence of the dipole-target amplitude is usually 
neglected, in order to disentangle the weak coupling dynamics associated with gluon saturation from the long-range 
strong coupling QCD dynamics. Going beyond that approximation remains the main challenge in gluon saturation 
theory |18] and will not be discussed further here. The imaginary part N{r,Y) of the dipole-target amplitude then 
depends only on Y and on the dipole size r. For physical solutions of the BK equation, N{r, Y) °^ in the dilute regime 
at small r and N{r,Y) — 1 in the dense regime at large r. At each Y, the transition between these two regimes is 
encoded by the saturation scale Qs{Y) and due to the BK equation, QsiY) grows with Y and N{r,Y) scales to a very 
good approximation with r^Qs^{Y) after a few units in Y. This is in qualitative agreement with the geometric scaling 
property found in the HERA data i 19il . 

The BK equation belongs to a class of nonlinear equations whose solutions have a universal asymptotic behavior, 
independent of the initial condition ll20ll2lll . In particular, the large Y behavior of Qs{Y) writes 

logQs^iY)=XY + ci logY + Const. + C2Y-^/^ + ff{Y-^) , (1) 

with A, ci and C2 three known universal constants 1221 l23i l24ll . However, for any realistic value of the coupling a^, 
the theoretical value A ~ A.E^NcCCs/n is much larger than the value A ~ 0.2 or 0.3 extracted from the HERA data. 
This too fast evolution of Qs{Y) despite the qualitative successes of gluon saturation shows the need for NLO/NLL 
corrections for a quantitative description of gluons saturation physics. 

RUNNING COUPLING AND OTHER NLO/NLL EFFECTS 

In order to reach the NLO/NLL accuracy, we need both the NLO contributions to the DIS cross section, and the NLO 
corrections to the BK equation so that it can resum high-energy next-to-leading logs. Both types of corrections have 
been calculated recently, in Refs. ll25ll and ll26l l27il respectively. However, more work is needed before a practical 
implementation of those results. Indeed, the NLO corrections to the DIS cross section are not yet available in the 
form generalizing the dipole factorization. And various large N LO corrections arising in the BK equation need to be 
resummed, in a similar way as for the BFKL equation 1281 l29l 130.1 . These large NLO corrections have three origins: 
the running of the coupling ttj, the kinematical approximations usually performed to calculate the high-energy LL and 
NLL contributions, and the non-eikonal parts of the DGLAP splitting functions. 

The large corrections of kinematical origin come from the incorrect ordering of successive gluons or dipoles in 
standard formalism. Leading logs indeed arise only when gluons are strictly ordered in fc+ and in simultaneously. 
Imposing this double ordering already in the LO equation then allows to resum large corrections of that type (sT]. A 
full treatment of these effects will be soon available [32]. Concerning the non-eikonal corrections to the BK equation, 
only a crude treatment is available |33t], so that more work is needed. 

The most studied and best understood of the large corrections are the ones related to the running of «(. A subset 
of the NLL contributions can be interpreted as LL contributions multiplied by the first correction in the perturbative 
expansion of the running coupling. These particular NLL contributions should not be included in the NLO corrections 
to the BK equation, but used instead to determine with what transverse scale runs in the BK equation. However, 
there is an ambiguity in this separation of NLL contributions in two types, which leads to various prescriptions for the 
running of the coupling |34, 35]. It is now widely accepted, both for theoretical and technical reasons, that Balitsky's 
prescription, proposed in Ref. 134.] , is the most appropriate. 

As an intermediate step between LO/LL and NLO/NLL accuracy, many studies have been dedicated to running 
coupling versions of the LO BK equation, in particular with Balitsky's prescription or more naive ones. In this context, 
a very good global fit of the HERA data has been performed 1.36.] . In that fit, however, Aqcd is taken as a fit parameter 
and is found significantly smaller than the physical value. 

Analytical studies ll22l |23|] have shown that when a, runs, the asymptotic behavior of the solutions of the BK 
equation is still universal, but in a different universality class than in the fixed coupling case. In particular, Qs{Y) now 



behaves as IHIIIIItI] 

\og{Qs\Y)/KlcD) = C,,2^ + Cy(,Y'l' + Co + C_i/6F-i/^ + C_,„Y-'I^ + ^(F-'/2) , (2) 

with Ci/2, Ci/6, Co, C_i/g and C_i/3 known and independent of the initial conditions. Moreover, C1/2, Ci/e C_i/6 
are totally insensitive to possible NLO and higher order corrections to the BK equation, whereas Co and C_i/3 depend 
on LO and NLO contributions, but not on NNLO or higher fyf\. Hence, at high enough Y , the main effect of NLO 
corrections to the BK equation apart from the running of is a shift of Co which can be absorbed by a rescaling of 
Aqcd- This explains the need for such a rescaling in the fit of Ref. ll36ll . where those NLO corrections are not included. 

The asymptotic behavior dU is consistent with numerical simulations of the BFKL equation with running as at 
LO, strict NLO, or NLO with resummations, with toy model gluon saturation effects included jssll . These numerical 
simulations also demonstrate the need to resum the large NLO corrections in order to obtain sensible NLO results, in 
particular at not too large Y . 
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